source('plotfdplnorm.R')

xs <- rep(0, 100)
ys <- rep(0, 100)

i <- 1
while(i <= 100) {
	xs[i] <- rlnorm(1, lmean,lsd)
	h <- M * lfdp(xs[i])
	ys[i] <- runif(1, 0,h)
	i <- i+1
}

points(xs, ys,
	pch='x', cex=.75,
	col='black', lwd=2)

